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1^ Abstract 

04 We study the geometric properties of a system initially in equilibrium at a critical point 

that is suddenly quenched to another critical point and subsequently evolves towards the 
f—{ new equilibrium state. We focus on the bidimensional Ising model and we use numerical 

fj methods to characterize the morphological and statistical properties of spin and Fortuin- 

^ Kasteleyn clusters during the critical evolution. The analysis of the dynamics of an out of 

C equilibrium interface is also performed. We show that the small scale properties, smaller than 

I the target critical growing length ^(i) ~ t^'^ with z the dynamic exponent, are characterized 

by equilibrium at the working critical point, while the large scale properties, larger than the 
critical growing length, are those of the initial critical point. These features are similar to 
what was found for sub-critical quenches. We argue that quenches between critical points 
could be amenable to a more detailed analytical description. 
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1 Introduction 

The relaxation dynamics of a macroscopic system taken to its critical point with some quenching 
protocol has received much attention. Most of the existing studies focused on the time evolution 
of global quantities (magnetization, susceptibility, correlation functions, etc.) and characterized 
their scaling properties with numerical simulations and renormalization group techniques [I] |2J 
[31I3]- These studies are not restricted to any spatial nor order parameter dimensionality. 

In two dimensional CTitical systems in equilibrium, powerful theoretical tools such as Coulomb 
gas methods [F, conformal field theory (CFT) [S^ and stochastic Loewner evolution (SLE) f7,'H",'5], 
allowed one to characterize the geometric and statistical properties of a large variety of mesoscopic 
observables in great detail. These objects give a more complete image of the system's equilibrium 
configurations than the global observables accessed with scaling arguments and renormalization 
group techniques. However, as far as we know, nothing is known about these objects during the 
out of equilibrium evolution of the same (and other) critical samples. 

An extensive numerical and analytic investigation of the coarsening sub-critical dynamics of 
two dimensional models from a mesoscopic point of view was carried out in recent years. The 
models treated were the clean Ising model with non-conserved |10|[TT1[T?| and conserved [T3] order 
parameter dynamics, the random ferromagnet with non-conserved order parameter dynamics [14j 
or still the q state Potts model [15, 16J . These studies allowed one to build a rather complete 
picture of the geometric and statistical properties of the spin clusters in these bidimensional 
systems. More precisely, their domain and hull-enclosed areas as well as their boundary lengths 
and the relation between areas and perimeters were analyzed and characterized in detail. 

The aim of this work it to present a similar study of relevant dynamic geometric objects 
during the critical non-equilibrium evolution of the 2d Ising model evolving from equilibrium 
at another critical point, in this case an infinite temperature configuration that is equivalent to 
critical uncorrelated site percolation. We use simple scaling arguments and extensive numerical 
simulations. 

Let us be more specific about the objects of our study. The most natural objects to consider 
are the spin clusters, i.e. clusters of nearest neighbor spins on the lattice that point in the 
same direction. These clusters are accessible via direct observation of the system. However, at 
the critical point, they are not appropriate to describe the critical equilibrium properties of the 
Ising model since their shape and statistical properties are not only dictated by the physical 
correlations but also by purely geometrical factors |171 [TFj . To remedy this problem one has 
to consider, in place, smaller clusters, namely the Fortuin-Kasteleyn (FK) ones, that capture 
exclusively the equilibrium physical correlations in the system. Although not relevant to describe 
the equilibrium physical macroscopic properties of the samples the spin clusters are, nonetheless, 
also critical at the phase transition. Consequently, they are characterized by a different set of 
exponents from the ones of the FK clusters. For this reason (to be explained more thoroughly in 



Sec. 2.2 1 spin clusters are also interesting to analyze and we study both types of clusters here. 

The article is organized as follows. In Sec. [2] we list the definitions of the different objects 
we study. The concrete analysis is presented next. We first check that the dynamical exponent 
z governing the critical dynamics of spin clusters is actually the one of the 2d Ising model, as 
obtained with other means, as the dynamic renormalization group method |2]. To the best of our 
knowledge, this has never been brought to light before. We then study in full extent the number 
densities of various quantities giving access to the structure of spin clusters on a large interval 
of sizes. The validity of the dynamical scaling hypothesis is tested upon these number densities. 
Relations between areas and boundaries are also explored. The same method is applied to the 
FK clusters. We later turn our attention to the dynamics of a single interface and the comparison 
with the equilibrium results given by conformal field theory and stochastic Loewner evolution. 
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This is the content of Sec. [3] Finally, in Sec. [4] we draw our conclusions and we discuss some 
lines of future research. 

2 Definitions 

2.1 Spin and Fortuin-Kasteleyn clusters 

Simple scaling arguments from percolation theory [19_ suggest that, in equilibrium, the divergence 
of the mean size of some kind of finite cluster at the critical point should be governed by the 
susceptibility exponent or, equivalently, the one of the spin-spin correlation function, that is 
the same as the probability for two spins to be in the same such cluster. The ensemble of 
nearest-neighbor spins that are parallel to each other constitute a domain or spin cluster. These 
are the most natural geometric objects, the critical properties of which one would expect to be 
linked to the ones of macroscopic physical observables. However, spin clusters do not capture the 
underlying critical properties of statistical physics models. This fact was first noticed by Sykes 
and Gaunt [17] who remarked that the mean size of the finite domains of the Ising model on a 
two dimensional lattice (2dIM) diverges at the critical point with an exponent that is different 
from the magnetic susceptibility exponent (see Table IlJ). Moreover, at infinite temperature there 
are spin clusters of arbitrary large size although the spin-spin correlation function vanishes. Still, 
for the particular case of the 2dIM, domains are critical at the critical point although they do 
not give access to the relevant critical exponents of the model. We will come back to this point 
in Sec. |2.2[ In three dimensions, this is not the case and these spin clusters do not even percolate 
at the critical point [20J. 

The point was then to build clusters containing just physical correlations between spins (and 
not trivial geometric ones). This is achieved by the Fortuin-Kasteleyn clusters [21', proposed by 
these authors as a way to describe with the random cluster model percolation, the Ising model, 
the Potts model and many other statistical models. Independently, Coniglio and Klein [TH] 
linked the critical properties of the Ising model to the percolation of Ising droplets which were 
later identified with the Fortuin-Kasteleyn clusters, the name that prevailed in the literature. 
The Fortuin-Kasteleyn clusters are constructed as follows. Starting with a spin domain, one 
first draws all bonds linking nearest-neighbor spin on the cluster and then erases bonds with 
a temperature dependent probability e~^^. In such a way, the original bond-cluster typically 
diminishes in size and may even get disconnected. More precisely, the construction of these 
clusters from the partition function Z is the following. One rewrites the exponential of the sum 
as a product of exponentials: 

with K = J/{kBT) — f3J, J the exchange coupling and ks the Boltzmann constant. Hereafter 
we measure temperature T in units of J/ks- By remarking that, since ai G { — 1; 1}, 

gif^.^, ^ ^K [(^ _ e~^^)S„^^,^ + e'^% (2) 

Z can be recast in the form: 

2=e"^"EnH...,+(i-rt], (3) 

with m the total number of bonds of the lattice (m — 3i~ for a planar triangular lattice of 
linear size L) and p = 1 — e^^^ . By expanding the product, if cr^ == aj one keeps the bond {i.,j) 
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with probability p or erases it with probabiHty I — p. The broken bonds (ai — —tj) are never 
kept. A connected collection of kept bonds is a Fortuin-Kasteleyn (FK) cluster. This whole 
construction is called the Fortuin-Kasteleyn random cluster model and an example is presented 
in Fig. [T] These clusters possess the searched properties. Their exponents are those of the Ising 
model [22 . For example, the probability for two spins to belong to the same FK cluster is 
exactly the spin-spin correlation of the Ising model. Note that the FK construction can be easily 
extended to all q-state Potts models, even with non integer q, in any dimensions. 






Figure 1: Sketch of the construction of a FK cluster. A domain is identified. The bonds between 
nearest-neighbor aligned spins (represented with black dots on the lattice sites) are erased using 
the FK procedure. Two FK disconnected clusters remain. The surviving bonds are highlighted 
on the edges of the lattice. 



2.2 A few hints on conformal field theory 

The existence of two types of critical clusters, the spin and FK ones belonging to distinct uni- 
versality classes is valid for all g-state Potts model with < g < 4 in two dimensions. In order to 
understand this point it is useful to consider the Coulomb gas formulation of the Potts model [5]. 
In this formulation, the 2d g-state Potts model can be described by the parameter k € [4, 8] such 
that 

y? = -2cos('^V (4) 

We choose a convention such that k also corresponds to the SLE parametrization for the interfaces 
associated to the g-state Potts model. We will come back to this point in Sec. |3.4| where we will 
study the behaviour of an out of equilibrium interface. The parameter n is also related to the 
central charge of the corresponding conformal field theory by the relation O E] : 



c = 



6- 



16 



(5) 



For K g [8/3,4] this Coulomb gas representation describes another class of models, the tri- 
critical Potts model, i.e. Potts models with dilution. Note that the central charge is invariant 
under the transformation k — >■ 16/ k which maps a critical g-state Potts model onto a tricritical 
Potts model with a different number of states given by eq. Q . This means that there exist two 
critical theories for a given central charge. One is associated to the critical Potts model and 
the relevant structures are the FK clusters. The other one is the tricritical Potts model and 
the relevant clusters are the domains or spin clusters. This duality is such that the geometrical 
clusters of one model are the FK clusters of the other model and vice versa |23l E^ . 

The tricritical model associated to the 2d Ising model {k — 16/3) is the dilute g = 1 Potts 
{k — 3). Both models possess the same central charge c = 1/2 but the same quantities in the 
two model are not associated to the same operators of the CFT. In 2d the spin clusters of the 
Ising model are the critical objects of the dilute g = 1 Potts model [55] . This explains why the 
percolation exponents associated to the spin clusters are not related to the Ising model exponents 
but in fact are those of the tricritical g = 1 Potts model, described by k = 3. 
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The critical exponents can then be expressed in term of k. For instance, the fractal dimensions 
Dc, Dfi and D^p of the cluster area, its hull and its external perimeter (a precise definition of 



these quantities will be given in Sec. 2.3| read |26|I23|: 

9 

(6) 

(7) 
(8) 

Note that the duality k — > 16/k relates the fractal dimension of the external perimeter to the 
dimension of the hull of the dual model, such as Depin) = Dh{16/K) for k > 4. This relation 
is called Duplantier duality [23]. For k < 4 the external perimeter coincides with the hull since 
they are no fjords. In this case the relation between Df.p and D^, mentioned above does not hold 
anymore and in fact D^p = D^. 



Dc -- 


^ -i-l 


Dh -- 


- 'n- 


Dep -- 


K 



Dc Dh Dep 



percolation 6 1 91/48 7/4 4/3 

Ising FK clusters 1/2 16/3 2 15/8 5/3 11/8 
Ising spin clusters 1/2 3 1 187/96 11/8 11/8 

Table 1: Central charge, c. Coulomb gas or SLE parameter, k, Potts model parameter, q [all 
these related by eqs. Q and (Is])], and three fractal dimensions at the percolation threshold, and 
for FK and spin clusters at the critical point of the 2rf Ising model. The c subscript in the fractal 
dimension is for the cluster mass, the h one for its hull, and ep for its external perimeter. 

Finally, let us note that the existence of two types of critical clusters at the critical point is 
specific to the 2d Potts model. While the FK clusters can be defined for more general models, 
nothing ensures that they are critical at the transition point. For instance, FK clusters are not 
critical at the critical point in 2d parafermionic models ^7}. Concerning the spin clusters, it is 
not even sure that they are critical at the critical point as in the 3d Ising model we have already 
mentioned. 

2.3 Definitions of the quantities computed 

In this section we define the quantities that we will consider in this work. We call a domain or 
spin cluster a connected set of sites with spins taking the same value. The mass (area in 2d) of 
a domain is the number of sites belonging to it. A broken bond is a link of the lattice between 
two neighboring spins with different value (see Fig. pi F] 

The domain wall of a spin cluster is its external and internal contour, constructed as follows. 
One first generates a dual lattice by placing a site at the center of each plaquette of the original 
lattice. Next, the links on the dual lattice that cross broken bonds on the original lattice are 
joined together. In this way, one finds a closed loop on the dual lattice that runs along the 
internal or external boundary of a spin cluster in the way sketched in Fig. [2] for the external 

^ In all the figures we use a square lattice for simplicity. The extension to the triangular lattice considered later 
should be straightforward. 
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# +1 Spin 

• -1 Spin 

— Broken bond 




— hull of length 



Figure 2: Left: sketch of an Ising spin configuration with the two- valued spins represented by 
black and red (black and grey) dots on the square lattice sites. Broken bonds are drawn with 
green (light grey) lines on the edges of the lattice. Right: a domain is singled out and its hull is 
represented with a thin red line. 



component. There is not much theoretical knowledge on this object and it is therefore more 
convenient to study other geometrical quantities which have received a theoretical description in 
equilibrium. 

The hull of a cluster is restricted to the external part of the contour, that is to say, one 
excludes the contribution of the holes of the cluster (see Fig. pi. The hull enclosed area is the 
area, i.e. the total number of sites, inside the hull (the holes within the domains are thus filled). 
In the example in the right of Fig. [2] the hull-enclosed area is 29. 

The length of the different types of contours defined above and living on the dual lattice is 
computed by counting the number of broken bonds corresponding to the object of concern that 
are crossed by the boundary. We can also define the external perimeter built by closing the 
narrow gates of the hull, making in this way a smoother version of the contour by eliminating 
the deep fjords. The meaning of narrow will be discussed later. As an example, on a square 
lattice, a cluster composed of a unique spin has a hull of length 4 and an area equal to 1 while 
for a two-spin cluster the hull length is 6 and the area equals 2. 

The FK clusters are defined on the edges of the lattice and not on the sites of the lattice as 
for the spin clusters. Therefore the length of their contour on the dual lattice is not proportional 
to the number of broken bonds. Indeed some bonds linking two sites of a FK cluster may not be 
within the cluster if they have not been kept in the construction of the cluster. In order to define 
the border of an FK cluster we then use a different sub-lattice, with four sites associated to an 
original one as shown in Fig. [3] The distance between two nearest-neighbor sub-sites counts as 
the unit of length for the contour. As an example, the cluster in Fig. [3] has a hull length of 24. 

In Fig. |3] we also present the construction of the external perimeter of an FK cluster. To 
measure its length, the procedure is the following: all the bonds between nearest-neighbor sites 
(on the original lattice) belonging to the FK cluster are drawn, and then the walker is allowed 
to wander around this new cluster, that is to say, on the sites of the sub-lattice introduced 
before. The external perimeter is a smoother version of the hull and it has, consequently, a 
smaller fractal dimension. The cluster in Fig. [3] has an external perimeter of length 20 which is 
effectively smaller than its hull length. 



2.4 Statistical and geometric properties in equilibrium 

To describe the statistical and geometrical properties of the spin and FK clusters we used the 
tools of percolation theory (see e.g. [19J. For example, we counted the number of spin clusters 
with a given area for numerous samples and we thus obtained, after normalization, the number 
density of the spin cluster areas. We followed this procedure for the hull lengths, the hull-enclosed 
areas, etc. In general, when we consider the probability distribution per lattice site n^lX) (of 
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perimeter 



Figure 3: Left: an example of an FK cluster. Center: the sites of the sub-lattice associated to 
the original square lattice are added as small points. Right: the hull and the external perimeter 
of the chosen FK cluster are shown with a thin green line. 



a given geometrical object x taking values X) which follows a power law, we denote the related 
critical exponent tx, i.e., 

n,{X)^X-^^. (9) 

Either the clusters or their boimdaries are fractal objects at the critical point. Their Hausdorff 
dimension (referred hereafter as fractal dimension) is non-trivial and has been related to other 
critical exponents. Indeed, the fractal dimension Dx of the objects described by the quantity X 
can be expressed in terms of tx ■ 

Dx = ^—. (10) 

Tx -I 

with d the dimensionality of the lattice [12] . A few examples are given in Table fll 

In finite-size lattices a cluster is said to percolate whenever it spans over a distance that 
is larger than the linear size L of the system, in at least one direction of the lattice. The hull- 
enclosed areas do not exhibit a fractal structure since they have no holes. Their fractal dimension 
is therefore 2. Their boundaries behave differently and they are fractal. Furthermore, Cardy and 
Ziff ^28j showed that the number density of hull-enclosed areas behaves as: 

nM) ^ ^ (11) 



A2 



for L^ > ^2 > a^ with C 



1/(8-^371") at the percolation threshold and C 
Ising spin clusters. The microscopic cut-off a is the lattice spacing. 



l/(16%/37r) for 



3 Critical dynamics 

The stochastic relaxation dynamics after a quench to the critical point has been performed with a 
time dependent extension of the renormalization group method. Scaling laws for averaged global 
observable such as correlation functions and others have been obtained in many different systems. 
The growth of an equilibration length, ^{t) ~ t^'^ with z the dynamic critical exponent, was 
evidenced. In this section we study the critical dynamics from a mesoscopic point of view and 
we show that this growing length also plays an important role in the characterization of the 
statistical properties of fluctuating quantities. 
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3.1 Description of the protocol 

In this section we study the evolution of a system after a quench from T = oo to Tc- We first 
consider spin clusters on a triangular lattice for which T^. = 4/ In 3. The choice of a triangular 
lattice is motivated by the fact that the infinite temperature point exactly corresponds to a site 
percolation critical point on this lattice. Indeed, the site percolation threshold for this lattice is 
Pc — 1/2 |29| so that, if we consider the number density of domains, the system is exactly at the 
percolation threshold at infinite temperature where the spins take the values 1 or —1 with equal 
probability. For the quench considered, the system is initially in equilibrium at a critical state 
and after the quench it evolves towards equilibrium at another critical state, the one at Tc- 

At infinite temperature FK clusters are not critical, actually they are trivial, for any lattice 
so it is immaterial to use the triangular or any other one. For this reason, in our study we used 
a square lattice. For the square lattice Tc = 2/ln(l + y/2) and the infinite temperature point is 
not even a critical one for the spin clusters since here pc > 1/2. (Having said this, the studies 
in [lUllTTlIT^ showed that, somehow surprisingly, the sub-critical Monte Carlo dynamics of such 
a non-critical initial state gets, after a few time steps, very close to the critical percolation point 
as far as the properties of spin clusters are concerned. For instance, their number densities 
rapidly develop the critical percolation tails and only later the dynamics evolve towards the 
target equilibrium state at low temperature.) 

The simulations were carried out on a lattice of linear size L with periodic boundary condi- 
tions in both directions. For the initial condition at infinite temperature, the spins were chosen 
randomly with equal probability of being up or down. Once the system was prepared in the 
desired initial condition, single spin updates were performed at the temperature of the quench. 
A Monte Carlo time step (MCs) corresponds to L^ single spin updates. The updates were ac- 
cepted or rejected via the standard Monte Carlo Metropolis scheme. We gathered around 6 10^ 
independent samples. Unless stated otherwise, the lattice used has a linear size L = 1000 for the 
analysis of spin clusters and L = 320 for the study of FK clusters. For each sample we computed 
the desired quantities every 2" MCs with n e [1, 12]. 

For the spin clusters, the algorithm distinguishes first between the internal and external part 
(hull) of the perimeter and then calculates the length of the hulls, i.e. the number of broken 
external bonds. To do so the lattice is scanned and when a broken bond that has not been counted 
yet is found, the algorithm follows the boundary in a precise direction and keeps track (with a 
cumulative angle) of the path followed. Then, depending on the sign of the angle the boundary 
drawn is an external or an internal one. Some clusters (spanning ones) have a boundary with a 
vanishing angle. Those boundaries run across the system from one side to another so they are 
not homotopic to a point on the torus. We checked that these clusters are sufficiently rare so 
that we can discard them without affecting the statistics. For the FK clusters the process is the 
same on the sub-lattice evoked above. The algorithm walks on the sub-lattice around the FK 
clusters and once a contour is formed, the angle is measured to discriminate between internal 
and external contours and their lengths are computed. 

The number densities obtained via this Monte Carlo method present a great dispersion, 
especially for large clusters which are much rarer than the small ones. This dispersion has been 
greatly suppressed by choosing appropriate bin sizes to construct the histograms, and to extract 
from them the distributions. 

3.2 Spin clusters 

We first briefly present the simulations performed to check whether the dynamical exponent 
governing the dynamics of the spin clusters conforms to the one extracted from the analysis of 
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the correlation functions that is well documented in the literature |30l [511 [3] . Then we present 
the results for the spins clusters on a triangular lattice after the quench. 



3.2.1 The dynamic exponent 



As explained in Sec. |2.2| the critical geometric clusters of the 2dIM are well understood as they 
are the FK clusters of the tricritical q = 1 Potts model. However, it is not obvious whether this 
correspondence should hold in an out of equilibrium situation, typically a quench towards the 
critical point. In particular, one may wonder whether the dynamical exponent z for spin clusters 
is the Ising one. To check this reasonable assumption, we calculated specific quantities that 
allow us to extract spin clusters exponents. For example, in equilibrium, the size of the largest 
spin cluster scales as L~(/^/'')t'i where {l3/v)tri — 5/96 and jitri is the magnetic exponent of the 
tricritical g = 1 Potts model, vtri is the exponent associated to the divergence of the correlation 
length and L is the linear size of the sample. We denote this quantity Mg. We now consider a 
quench from T — Q to Tc for different system linear sizes L, and we compute Mg{t). We consider 
this quench since the scaling form is simple in this case and the dynamical exponent should not 
depend upon the initial condition as long as we quench at the critical temperature. The inset of 
Fig. HI shows the relaxation of this quantity for several sizes. We then apply the following scaling 
Mg{t) — )■ Mg{t)L'^^/^^*" and t — >■ t/L^ . The collapse of the different curves giving access to z is 
presented in Fig. |4] The value obtained converges towards the most accurately estimated value 
z = 2.1667 (5) [31j. This gives strong evidence that the dynamics of the spin clusters are indeed 
governed by the Ising dynamical exponent z. In consequence we will use the same z for all the 
quantities computed in this work. 




t/L" 



Figure 4: Geometrical cluster magnetization Mg{t). The main part shows the collapse onto a 
master curve with the rescaling t — > t/U^ and Mg{t) — > AIg{t)L^-^/'^^^^\ In the inset we show the 
raw data Mg{t) versus time for different system linear sizes L. 
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3.2.2 Hull-enclosed area number density 

In Fig. Is] we present the distribution of the hull-enclosed areas, nh{A, t), at different times after 
the quench. In the first panel. Fig. p[a), we display the distributions, ri/j, versus the hull- 
enclosed area A. As expected from the results of Cardy and Ziff [2S], the slopes of the initial 
and asymptotical (equilibrium at Tc) distributions equal —2 in a double logarithmic plot but 
the prefactors are slightly different. For this reason, the changes occurring in the distributions 
are very small. Equilibrium data (generated with a special purpose algorithm at Tc) are shown 
with (green at T — >■ cx) and red at T^) dashed lines. The bumps at large sizes correspond to the 
spanning clusters with a linear size of the order of L, the linear size of the system, are finite size 
effects. Note that the position of the bump and the very last part of the n/,,(A, t) close to it, are 
time-independent (within our numerical accuracy and for the times accessed in the simulation). 
For hull-enclosed areas dynamic scaling suggests 

n^{A,t)^A--g(^^^ (12) 

for Aq -^ A -^ L^ with Aq — a^ a, microscopic area scale and L^ the macroscopic one. To make 
the notation lighter we did not add any sub-script to t and a. The exponents are given by 

T = l + d/D, a:^D/z. (13) 

Since these areas are regular, i.e. they have no holes, D = d = 2 and 

r = 2, a==2/z. (14) 

Consistency with the asymptotic limit requires g{y -^ 0) = C and g{y — > oo) = 2C. 

In order to test the dynamical scaling hypothesis, we impose the value t — 2 and rescale the 
areas by a factor t" with a a free parameter whose value is determined by the collapse of the 
numerical data on a master curve. The best collapse is obtained for a — 0.92 (5) and is seen 
in Fig. plb) with the theoretical percolation and Ising critical distributions being just horizontal 
after the given ordinate rescaling. The numerical equilibrium critical distributions (red dashed 
line for Tc and green dashed line for T — )■ c») have been placed at convenience (horizontally) 
to ease the comparison with the other distributions. The horizontal dotted lines correspond to 
the exact prefactors C = l/{lQ^/3n) and 2C for the critical Ising and percolation distributions, 
respectively. The value of a is in agreement with the expected value a — 2/z ~ 0.92306 with 
z ~ 2.1667 (5) |3ij since the hull-enclosed areas scale as the square of the dynamical length scale 
S,{t) ^ t^l^ . It is quite remarkable that, while no time-dependent rescaling has been applied to 
the distribution, the curves collapse so accurately in the vertical direction. One may remark that 
the scaling is not so good for times that are smaller than 64 MCs, which is not surprising since 
the dynamical scaling hypothesis only holds after a non-universal time-scale. 

In curvature-driven coarsening we know from [TTl [T^ . that n/j(A,t) = 2C/(yl -I- A/jt)^ if we 
use an infinite temperature initial condition, i.e. n/j(A, 0) = 2C/A^, with A/j a non-universal 
parameter. It implies 

on 

A^nh{A,t)^f{A/\Ht) with /(a;) = ^^-^^. (15) 

f{x) looks like the right part of A'^nh{A, t) for A > t^/^ in our case. Since we consider a situation 
with two critical points and not only one as in the sub-critical case we are tempted to think that 
the minimum of g{x) around a; ~ 1 and more generally its non-monotonic behaviour can be 
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Figure 5: Spin clusters hull-enclosed area number densities in the 2dlM on a triangular lattice 
are presented, (a) With dashed (red and green) lines equilibrium data at T — ;> oo and T^- With 
solid lines the raw dynamic data for the times given in the key (in MCs) including the spanning 
clusters using a double logarithmic scale, (b) The areas in the x-axis are rescaled by the factor 
t", with a = 0.92 consistently with the expectation a = D/z = 2/z. The distributions in the 
y-axis (in linear scale) are multiplied by A'^ = A^. The grey dotted horizontal lines at 0.046 and 
0.023 are the theoretical predictions for the T — >■ oo critical site percolation point and the T^ 



one, respectively, and g(x) is the fitting function defined in eq. (12 1, and discussed in the text 
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reproduced with a sum of two functions similar to f{x), one decreasing from C to and one 



increasing from to 2C. This suggests for g(x) defined in eq. (12) the form 



gi^) = C -, T— + -, -, T—-r^ . (16) 

We left the powers h and c as free parameters as there is a priori no reason for them to be simple 
integers as in the sub-critical coarsening situation. Numerical inspection of the data shows a 
good fit for the values a ~ 0.65, h ~ 0.79 and c ~ 2.3. Roughly speaking, a fixes the position of 
the minimum, h its width and c its depth. Note that c is chosen to allow the fitting curve g(x) 
to go through the minimum of the numerical data but we do not attribute a special meaning to 
this value. Indeed the expansion of g(a;) close to the asymptotic values C and 2C is proportional 
to x^ and x~^ , respectively, both independant of c. 

3.2.3 Hull-enclosed areas and hull lengths 

As in pn [T3J [15] we study the relation between the hull-enclosed areas and the hull lengths by 
tracing an averaged scattered plot of A against p in Fig. p] In panel (a) we show the raw data 
for different times after the critical quench. In panel (b) we scale the data by the relevant typical 
growing scales. In the case of the hull-enclosed areas this is the linear growing length, ^(i), to 
the power of their fractal dimension that is simply D — d=2iox these regular objects. In the 
case of the hull lengths, instead, the relevant linear scale is still the growing length, ^(i), now to 
the power of the hulls fractal dimension in the critical Ising point, which is D\^ = 11/8, as given 
in Table [T] The master curve shows a clear cross-over between two power law behaviors with the 
powers D / Dj^ = 16/11 ~ 1.45 controlling the small scales and D/Dj^ ~ 8/7 ~ 1.14 controlling 
the large scales. Summarizing, 

A_\ ^ f^L_V ^ith C = I ^{^i ''' '"T, « ^' (17) 



tD/zJ \tD'J'-J ^ \ D/Dj^ for p/i^'^/^>l. 

3.2.4 Hull length number density 

Next we proceed with the hull length distribution displayed in Fig.lTJ In the first graph, Fig.l7[a), 
we show the distributions nh{p,t) vs. p for various times after the quench where p stands 
for the hull length. For this quantity, the equilibrium behavior for the critical Ising model is 
"i^hip) ^ p^27/ii g^^^ j|. jg different from percolation criticality where nh{p) ^ p^^^l'^ . These 
equilibrium curves are drawn with dashed lines in the figure. We see in Fig. JT^a) that the 
dynamic curves interpolate between these two critical laws. For example, we can observe that for 
t = 4096 MCs there is a qualitative change for p ~ 2000. However, since the difference between 
27/11 and 15/7 is small it is still difficult to get a precise picture of what is happening with this 
data representation. Note that, contrary to what happens with the area number densities, the 
time-dependence is seen in the full extent of the curves and not only for (relatively) small scales. 
Even the bump is displaced towards shorter lengths in the course of time. 

A better description is given in Fig. ^b) where we drew the same distribution multiplied by 
p^ with T = 27/11, the exponent of the tail in the equilibrium distribution at Tc, so that the 
Ising equilibrium critical distribution is horizontal. As before, the lengths are rescaled by a factor 
i" to take into account the growing length scale. Searching the value of a that gives the best 
horizontal collapse we find a = 0.63 (3). This fits well with the value expected from dynamical 
scaling argument such that p is rescaled by ^{t)^'^ with D/, = Dj^ = 11/8 = 1.375 the fractal 
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Figure 6: The hull-enclosed area against the length of the hull at different times after the critical 
quench for L ~ 400. (a) Raw data, (b) Scaled plot. The master curve crosses over from x'^'^'^ 
to x"^/^^ with d/Dl ~ 1.45 and d/Df^ ~ 1.15. 
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dimension of Ising clusters' hulls. D^/z ~ 0.6348 is in good agreement with the measured 
value and this clearly supports the idea that the growing clusters have a fractal boundary. It 
is interesting to notice that this does not happen in sub-critical coarsening where the domain 
growth is curvature-driven and the boundaries of clusters that are smaller than the cross-over 
scale are smooth and do not have a fractal structure (see [TT] e.g. for the study of such a case 
in the 2dlM). The asymptotic value, for p/f-^'^^^ <C 1, coincides with the expected equilibrium 
value, as obtained from the equilibrium simulations and shown with the (red) dashed horizontal 
line. 

The behavior in the large scale or short time limits, p ^ t^^/^^ is also interesting. Since 
there is a one-to-one relation between areas and perimeters, we can use nfi{A,t)dA ~ nfi(p,t)dp 
where A and p are related by eq. ( [T7| (note that we use the same symbol Uh to represent two 
different functions). Using nh{A, t) ~ A~'^ with t = 1 + d/D for large values of A and the value 



of the exponent C, given in eq. (17 1 for large scales we find 



nu{p,t)p< ^ [-^y with j^d^^p^=rl-r^ (18) 

independently of D. We used the super-script / to stress that the r- value in the factor in the left- 
hand-side corresponding to critical Ising equilibrium was used in the scaling of data in Fig. ^b). 
The bump at the end of the data is just the contribution of the spanning clusters. The growing 
part of the master curve is very well described by the power law given above. Indeed, close 
to it we placed the data obtained in equilibrium at T — > c» (dashed green line) that in the 
representation used in the plot is given by x'', with the value of 7 given in eq. ( |18[ ). The two 
curves are parallel in the selected range of variation (within numerical accuracy) confirming our 
prediction. 

3.2.5 Domain area number density 



We have checked that the area- hull relation for the domains conforms to the scaling in eq. (17) 
with the relevant C given by the fractal dimension of the domains. 

Finally, we present the spin clusters (or domain) area distributions at several instants after 
the quench in Fig. ^ Again the two equilibrium distributions appear as (green and red) dashed 
lines. Even for the relatively short time t = 8 MCs the change in the distribution is rather 
drastic, and we can observe that the system has been depleted from small clusters with area 
up to 80 because of the onset of the interactions. This is related to the fact that the fractal 
dimension of the Ising clusters at equilibrium is greater than the dimension of the percolation 
clusters. Qualitatively we can explain this fact by saying that the Ising clusters are less porous 
because of the interactions. We perform again the rescaling of the areas by t" and we multiply 
the distributions by A'^ , with r = 379/187 ~ 2.027 being the critical Ising exponent for the 
distribution of domain areas (see the Table). Here the best collapse (Fig. Islb)) is obtained for 
the value a = 0.90 (3). The fractal dimension of the geometrical clusters at the critical point 
is Dc = 187/96 ~ 1.948. Then the value of a is in good agreement with the value expected 
from scaling arguments Dc/z ~ 0.899. Note that if we had neglected the fractal character and 
used Dc = 2, we would have obtained 2/z ~ 0.923 which is also in the interval of confidence 
of our numerical estimate but we think the value Dc/z with D — 187/96 is the correct one. 
The scaling just discussed gives evidence for the importance of the fractal structure of the spin 
clusters at the critical point. The values obtained suggest a fractal front, whose fractality is the 
critical Ising one, propagating on increasingly larger scales and bringing the system towards a 
new equilibrium, although very slowly, with a law governed by the dynamic exponent z. 



3 CRITICAL DYNAMICS 



16 



Ising 
percolation 



' — n- 

t=cx> 

t=0 
8 

64 
512 
4096 





Figure 7: Time-dependent distribution of the hull length p of the spin clusters, in a double 
logarithmic scale for the times given in the key in MCs. (a) Raw data dynamic data together 
with the ones for equilibrium at the departure and arrival values, T — ^ oo and T^, that are shown 
with dashed (red and green) lines, (b) Same as in (a) but the distributions are multiplied by 
p^ with T = 27/11 the exponent of the equilibrium number density at Tc and the lengths are 
divided by t" with a = D^/z ~ 0.63 with D^ the fractal dimension of the equilibrium hulls at 

Tr. 



3 CRITICAL DYNAMICS 17 



3.3 Fortuin-Kasteleyn clusters 

In this Section we present the results of simulations in which we analyze the FK clusters. These 
clusters are critical at T = T^, but they are not critical at T — )■ oo since at infinite temperature 
p — 1 — e~^^ = and they are all composed of only one site. We constructed the dynamic 
FK clusters adapting the usual procedure explained in Sec. |2.3| and sketched in Fig. [T] to the 
time-dependent domains. At different instants after the quench from /3 = to /3c we used the 
probability p — 1 — e^'^^'= to erase bonds from the spin domains existing at the chosen times. We 
next studied various distribution functions associated to the FK clusters. We recall that in this 
study we used a square lattice for simplicity and that the initial point (infinite temperature) is 
not critical for the spin clusters either. 

We have checked that the number density of hull-enclosed and domain FK areas as well as the 
area-perimeter relations satisfy the expected scaling relations already discussed for spin clusters 
in Sec. |3.2| so we prefer to show the data for other observables. 



3.3.1 Hull length 



In Fig. 10 a) we show the hull length distributions at T — > oo and at Tc as well as the distributions 
of the same quantities at several times after the quench. We can easily see that the initial 
distribution is not a critical one since it is not a straight line in this log- log graph. As time 
elapses, the distribution approaches a power law behavior since it heads towards criticality. 
Again, to get a better description, the curves are collapsed on Fig. [To|b) by using the now usual 
rescaling p^nfi{p,t) = /{p/t"). The best collapse onto a master curve is obtained for 6 = 2.12 
and a = 0.87. These values are compatible with r = 2.2 and D^/z ~ 0.77, the values at the 
Ising critical point, cf. Table [l] though here the agreement is not as good as for the hulls of the 
spin clusters. We notice that for times longer than 2000 MCs the curves do not collapse anymore 
for large sizes. This is expected since it roughly corresponds to the moment when the dynamical 
length scale reaches the order of the linear size of the system ^(i) ^ L. After this time we are no 
longer in the regime of unconstrained domain growth as the clusters are sensitive to the finite 
size of the system. 

3.3.2 External perimeter 

The same scaling is finally applied to the external perimeter distributions of the FK clusters, the 
idea being to check whether the consequences of Duplantier's duality hold out of equilibrium. 
The curves are not shown since they are very similar to the ones in Fig. [TO] The curves collapse 
for 6 = 2.56 and a = 0.63. To compare this values with the related equilibrium exponent at Tc 
we must remember that, according to the duality relation evoked before [Dep{K) = Dh{lQ/K) 
for K > 4] , the external perimeters of the FK clusters at T^ should scale at equilibrium like the 
hulls of the spin clusters, that is possess also at Tc a fractal dimension Dcp{K = 16/3) = Dh{K = 
3) = 11/8. This gives r ~ 2.45 which is compatible with S and a = Df,p/z ~ 0.63. Therefore the 
relation between the statistical and geometric properties of the external perimeters of FK clusters 
and the hulls of spin clusters proven in equilibrium, seems to remain valid out of equilibrium for 
the growing FK clusters. 

3.4 The interface 

In the previous sections, we have seen fractal dimensions appear while considering the dynamic 
scaling of several probability distributions. A different strategy, that has proven to be very useful 
in the analysis of the equilibrium critical points and in particular in the context of the study 
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Figure 8: Time-dependent domain area distribution, (a) Raw data at the times given in the 
key in MCs. (b) Same as in (a) with the distributions multipHed by A'^ — yl2.025 ^^^^^ ^j^g ^^^g, 
rescaled by t" with a = D/z ~ 0.90. 
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of SLE [71 [SI [3, is to study the fractal properties of an artificially generated interface. This is 
the line of research we follow here by generating the interface as follows. We take a 2dIM on a 
triangular lattice and we quench the system from equilibrium at T —?► cx) to Tc at i = 0. With 
appropriate boundary conditions we force the existence of a unique curve, defined on the edges 
of the dual honeycomb lattice, separating spin clusters of opposite sign and going, say, from top 
to bottom. One such curve is shown in Fig. [9] their properties are analyzed in this section. 




Figure 9: An interface in a 2dIM on a triangular lattice created by imposing spin up boundary 
condition on the left (grey hexagons) and spin down on the right (white hexagons). The curve 
is defined on the edges of the dual lattice, the honeycomb lattice represented in the figure. 



3.4.1 Fractal dimension 

We are interested in the fractal dimension of the interface. To compute this quantity we measured 
the length of the curve 1{L, t) with L the linear size of the system as a function of the time t. If 
the curve is fractal at time t, then the fractal dimension df{t) is defined by l{L,t) ex L'^-'^*). An 
effective df{t) is obtained from the slope of l{L,t) between two sizes L and L': 

^/(^'^'')- InL-lnL' ' ^^^^ 

This yields a fractal dimension depending on time but also on the two lengths L and L' because 
of finite-size corrections. In Fig.fTTlwe draw this fractal dimension vs. time t rescaled by L^ for 



several values of L and L' with z the dynamical exponent employed in Sec. 3.2 and Sec. 3.3 

We observe that for short times, i.e. t <^ L^ , the fractal dimension is compatible with 
7/4 = 1.75 the fractal dimension of the hulls at the percolation threshold, while for long times it 
is closer to 11/8 = 1.375 the fractal dimension of the Ising clusters' hull at the critical point. The 
compatibility is strong since even the finite-size corrections for the short (long) scales correspond 
to the related equilibrium values. The points completely on the left correspond to the initial 
values being rejected to — oo in the log scale. The dynamical scaling is also very accurate here 
since the curves perfectly collapse apart from finite-size corrections. 

The main disadvantage of calculating a fractal dimension this way is having to handle the 
two lengths in the interpretation of the fractal dimension. This is uneasy. We use a different and 
more performant approach in the next subsection. 
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Figure 10: Number density of FK hull length p at several times after the quench given in the 
key for a system of linear size L = 320. (a) Raw data, (b) Same data presented in the rescaled 
form p'^nf^^ {p,t) vs. p/t" with t = 2.12 and a — 0.87 that gives the best data collapse. 
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Figure 11: Fractal dimension of an interface vs. t/L^ 
of system sizes L and L' — 2L given in the key. 



as deduced from eq. ( 19 1 for several couples 



3.4.2 The winding angles 

The percolation and Ising model hulls, at equilibrium and in the continuum limit, are conformally 
invariant curves described by a stochastic Loewner evolution SLE^ with k = 6 for percolation [35] 
while K = 3 for Ising hulls [331 [34j- The parameter k can be determined numerically by computing 
the variance of the winding angle for two points chosen at random on a conformally invariant 
curve at a distance I along the curve which equals 



= ct 



4P - 1) 
D 



\nl, 



(20) 



where D is the fractal dimension of the curve [551 1551 [57] and ct is a constant. Finally, using 
eqs. (It]) and ( 20 1 , we obtain the following form for the variance of the winding angle: 



= ct 



4k 



InL 



(21) 



The winding angle variance, {6'^){l,t), is measured as a function of the distance / for different 
times after the quench. The results for a system of linear size L — 1280 are given in Fig. fl2ja). 
For the time t — Q corresponding to percolation, a fit of the form (21 1 yields k 
agreement with the exact result k, 



5.98 in excellent 
4096 MCs. a fit of the 



form (21 1 yields k 



6. For the largest time simulated t 
2.96. For this fit we kept the data with h\l < 6.5 since for larger distances, 
the variance starts deviating from a straight line. This value is in excellent agreement with the 
exact result for Ising model spins clusters in equilibrium, k = 3. The best fits for i = MCs 
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Figure 12: Variance of winding angle of interfaces out of equilibrium (a) and in equilibrium (b) 
for a system of linear size L = 1280. (a) (6'^) vs. In? for different times after the quench from 
T = oo to T = Tc- (b) (0^) vs. In? at equilibrium for different inverse temperatures /3 — \/T in 
unit of j3c given in the key. 
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and t — 4096 MCs are represented as dashed lines in Fig. [I2[a). We now wish to collapse these 
curves to make the dependence of the growing length scale upon time explicit. The length I 
along the curve is divided by t" and the winding angle variance corresponding to a length t" 
is subtracted. Thus, in Fig. [is] we plot {e'^){l,t) - 4k/(8 + «;)lnt" vs. ln(t/i")- The curves 
collapse very well for a = 0.65 which is in good agreement with D^/ z ~ 0.63 with D^ — 11/8 
the fractal dimension of Ising clusters' hulls. This means that the same rescaling as the one used 
for the distributions works for the winding angle variance as well. In the latter case, however, 
the transition between the two equilibrium regimes is more rapid. Indeed, for distances smaller 
than the characteristic length ^(t)'^'' ~ t'^h/z ^.j^g interface is, at least for the observable (6'^), as 
if the system were in equilibrium at T^ and for larger distances as if the system were still in the 
independent percolation regime. In the previous section, we have observed a similar crossover in 
the behavior of the probability distributions, with small and large scales being the ones of the 
initial regime of percolation and the ones of the state towards which the system is converging, 
respectively, but the crossover was far less abrupt. 

For comparison, we present the plot of {O"^) against In^ in equilibrium for various temperatures 



in Fig. 12 b). It is known that the high temperature phase of the Ising model at zero field on 
the triangular lattice is peculiar in the sense that it is critical and in the universality class 
of uncorrelated percolation [551 [55] . This is not true for other planar lattices. As we have 
already seen, the contours of percolation clusters on the triangular lattice are described by 
SLEg. This explains why for sufficiently large sizes (larger than the correlation length induced 
by the interactions between spins) we recover a dependence of (0^) with InZ that is compatible 
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with K — 6, i.e. the behavior of the interface is governed by the attractive infinite temperature 
point. 

Both Figs. 12 'a) and [T2[b) are compatible with k = 6 at large scale. However, for shorter 



distances the behavior of the quenched and equilibrium curves are very different. In the case of 
the quench, for sufficiently small size and up to Z ^ f^h/z ^.j^g curves are compatible with k = 3 
(Ising spin clusters) but for the curves at equilibrium the critical Ising behavior is not observed 
at all. 

Thus during the coarsening process the interface has a radically different behaviour from the 
one obtained by considering data for equilibrium at indermediate temperatures. In the first case, 
we flow from T = oo to T^, with an increasing scaling length ^(t) ^ t^'^. In the second we flow 
from Tc to r = oo with a correlation length £, ^ \T — T'cl "^ with 1^ = 1. In this second case 
we can only observe the critical large scale behaviour dominated by the percolation point while 
in the first case we distinguish two behaviors with a crossover controlled by the characteristic 
length $(t) - t^/\ 

4 Conclusion 

In this paper we analyzed the dynamics of a system instantaneously quenched from one critical 
point to another one. We used the prototype statistical model, the 2d Ising model, quenched 
from T — > cxD, i.e. critical site percolation, to the Ising critical point. 

As observed in [in]-[TB] for sub-critical quenches, the typical growing length, ^(i) ~ t^/^ in 
this case, separates two length scales. In the smaller one, all mesoscopic observables (areas, 
perimeters, etc.) satisfy the statistical and geometric properties of the working temperature 
equilibrium state. Instead, in the largest scale, these same observables are the ones of the initial 
state. In sub-critical quenches, this result could be shown analytically for the hull-enclosed area 
distribution [101 113] and it was confirmed numerically for different quantities and in different 
systems. 

The attraction of quenches between critical points is that, in the departing and arrival equi- 
librium conditions, very powerful analytic methods (Coulomb gas, conformal field theory, SLE) 
allowed one to compute a myriad of geometric properties including fractal dimensions and re- 
lated critical exponents. Could these methods be extended to deduce, analytically, the results 
presented in this manuscript for a particular case, the 2d Ising model, that we conjecture apply to 
all critical quenches, at least in bidimensional systems? This is definitely an interesting question 
on which we plan to work in the future. 
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